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^ I The next supernova in our galaxy will be detected by a variety of neutrino detectors. In 

O ' this lecture I discuss the set of observables needed to constrain the models of supernova neutrino 

I emission. They are the flux normalizations, and average energies, of each of the three expected 

components of the neutrino flux: Vf., S'e, and Vx (all the other four flavors combined). I show how 
the existing, or soon to be operational, neutrino detectors will be able to determine the magnitude 
of these observables, and estimate the corresponding rates. 

> 

1 Introduction 

When neutrinos from supernova 1987A were detected by the Kamiokande Q and 1MB Q collabora- 
tions, a new era of neutrino astrophysics began. Despite the limited statistics (11 events in Kamiokande 
and 8 events in 1MB) the observation confirmed that the core collapse supernovae emit most of their 
binding energy (a few x 10^^ erg) in neutrinos, that the duration of the neutrino emission is ~ 10 
seconds, and that the average energy of the neutrinos (at least of the v^, which were the only fiavor 
O \ actually seen) is ~ 15 MeV, close to expectations. The observation of SN1987A lead to a flood of 
^ \ papers analyzing its consequences (for a relatively early review, see e.g. [^), which is only very slowly 
diminishing with time. 

Historically, there were seven supernovae in our galaxy proper recorded in the past thousand years, 
^ \ and none in the last three centuries (some were not core-collapse SN which emit neutrinos, though). 
' All of them were relatively close to the solar system, so it is difficult to estimate the true rate averaged 
over the whole galaxy from this record. Consensus estimate of core-collapse supernova rate in our 
galaxy is about three times per century Q. Thus, the next Galactic supernova neutrino burst can 
come at any time, tomorrow or in several decades. It is likely that neutrinos from such supernova will 
be detected by a variety of detectors, with much better statistics than for SN1987A. Thus a wealth 
of new information is expected from such unique event which cannot be repeated in the productive 
lifetime of an average physicist. (Unfortunately, the present or planned neutrino detectors are unable 
to observe supernovae in even the nearest galaxy, Andromeda, about 700 kpc away.) Here 1 discuss 
some of the lessons that should, and hopefully will, be extracted from the neutrino signal of the next 
supernova in our galaxy. 

There are several areas of physics that will greatly benefit from the supernova neutrino observations. 
They can be divided into three broad categories: 

1. ) Neutrino properties; mass, mixing, decay, etc. In particular, one could use the time-of-fiight of 
the neutral current signal (dominantly and v^-) to reach sensitivity to masses of about 30 eV 



*Talk at the School 'Neutrinos in Astro, Particle and Nuclear Physics', Erice, September 18-26, 2001. 



1 



for these neutrinos This would represent an improvement by more than three orders of 

magnitude for the mass associated with f^, and by almost six orders of magnitude for the mass 
associated with when compared to the present direct neutrino mass limits [0]. If, moreover, 
the neutrino emission is abruptly truncated by the collapse of the proto-neutron star into a black 
hole, one can use this sharp cut-off in the neutrino signal to improve the time-of-flight sensitivity 
to masses of ~ 6 eV for and Vj- and to ~1.8 eV for v^. 

2. ) Supernova properties. From the neutrino signal it might be possible to determine the lumi- 
nosities and average energies of all three components of the neutrino flux: v^,, D^,, and z^^- (this 
notation will be used from now on collectively for i/^, v^- and their antiparticles). 

3. ) Supernova localization. Using the angular distribution of the products of the neutrino induced 
signal, or the timing of the signal recorded in widely separated detectors, it might be possible 
to find the direction towards the supernova independently, or prior to, of the optical signal (for 
the discussion of this item, see 1^]). 

I refer to the listed references regarding the items 1.) and 3.) and in the following I will concentrate 
on the item 2.) - the determination of supernova properties from the neutrino signal. My aim is 
going to be a definition of a 'template', i.e., a recipe how to determine the required quantities and 
what signal and statistical accuracy one may expect using the existing, or soon to be operational, 
detectors. Substantial deviations from this template will mean either that the supernova behaves 
in an unexpected way, or that neutrino oscillations affect the signal. Obviously, general analysis of 
all possibilities is impossible before the fact. However, the existence of such a template might help 
in preparing the detectors for the supernova signal, particularly those like SuperKamiokande, SNO, 
KamLAND or Borexino, which are built for a different purpose. 

I will consider a 'standard' supernova (for a review of Type-II supernova theory see approx- 
imately at the center of the galaxy, at the distance from Earth of 10 kpc. The binding energy, which 
is essentially fully emitted in neutrinos, is assumed to be 3 x 10^^ ergs. It is easy to understand the 
magnitude of the binding energy Eb^Y using the simple estimate 

Eb ^ , where R = IQ km, M ~ IAMq . (1) 

5 R 

Neutrinos are trapped in the hot and dense protoneutron star. The mean free path of neutrinos, 

A = — ~ 10 m for p ~ 10^*^ nucleons/cm^ , a ~ 10"^^ cm^ (2) 
pa 

is substantially shorter than the radius of the protoneutron star. In fact, the trapping occurs already 
when the star radius is ~ 100 km and the mean free path becomes comparable to the scale height h 
of the infalling matter {h = kT/MpQ where g ~ 10^^ ms~^ is the gravitational acceleration at that 
radius). 

Trapped neutrinos diffuse through the protoneutron star. They leave the star when they reach 
the so-called neutrinosphere, essentially the radius where their mean free path is comparable to the 
corresponding scale height at that point. Again, a crude estimate of the diffusion time is just the 
product of the time duration between successive scatterings and the number of steps, 

A 

^diff ~ ~ 10 s . (3) 

Throughout the star, neutrinos of all flavors are in equilibrium, with decreasing temperature at 
increasing radii. Thus, the temperature of the outgoing neutrino flux for each flavor will be the char- 
acteristic temperature of the corresponding neutrinosphere. Since the mean free paths of the different 
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neutrino flavors are different, the position of their neutrinospheres, and hence also the decoupling 
temperatures, will be different as well. The Vx neutrinos undergo only neutral current interactions, 
hence their mean free path is longest, and thus their decoupling temperature will be highest. Both v^, 
and v'e have in addition also charged current interactions. Moreover, since the star contains many more 
neutrons than protons, the mean free path will be shorter (since z/g interact with neutrons) than 
the v'e mean free path (since v'e interact with protons). Hence a hierarchy of decoupling temperatures 
(or mean energies) is expected, 

r(zy^)(~ 8MeV) > r(z^e)(~ 5MeV) > T(i/e)(~ 3.5MeV) , (4) 



or 

{E^J ~ 25 MeV > {EpJ ~ 16 MeV > {E^J ~ 11 MeV . (5) 

At the same time, one expects that the total luminosity will be equally shared by all neutrino 
flavors, so averaged over time 

Eb 

{Ly) ~ ~ 5 X 10^^ erg/s for all 6 flavors . (6) 

Note that the the short initial v^, neutronization pulse has only small luminosity when compared to 
{Ly) and is going to be difficult to observe. For a detailed description of the supernova neutrino 



emission, including the justification of the choice of the decoupling temperatures, see Refs. |11, 12 1. 
For each neutrino flavor the corresponding time averaged fiux at Earth will be therefore 

2.6 X IQll _2 -1 

^^-^ = {Ey){MeY) ^ ' 

for the assumed 10 seconds emission time. With such a flux and a typical cross section of ~ 10^^^ 
cm^, one expects few hundred charged current interactions with protons in 1 kton of water, and few 
tens of events in 1 kton of iron (or other heavy target). Clearly, very large detectors, operating for a 
long time, are needed. 

Thus the challenge for supernova neutrino observers is to detect separately the three expected 
components of the neutrino flux: the v^, component through the charged reaction on bound neutrons 
(i.e., on nuclei), the component most easily through the charged current reaction on free protons, 
and the Vx component through neutral current reactions. For each of these components one should 
determine, ideally, not only the total rate, proportional to LyjlEy) J a{Eu)f{Eu)dEu, [f[Ey) is the 
normalized energy distribution; typically assumed to be the Fermi-Dirac thermal one) but also the 
temperature, or equivalently {E^). If, and only if this program can be accomplished, can one reach 
reliable conclusions about supernova astrophysics and/or neutrino oscillations. 



2 Detecting Up and z/g through charged current reactions 

It is relatively easy to detect Ve-, since most detectors contain free protons and one can utilize the 
reaction + p e"*" + n with large cross section and a characteristic signature of the time and 
position correlated positron and neutron. 

The cross section is well known. Neglecting the small neutron recoil energy (~ E'^/Mp), one can 
simply relate the positron energy to the incoming neutrino energy, 

=Ey-A, A = Mn-Mp = 1.293 MeV . (8) 
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The differential cross section to this {Mp — > oo) order is 

^ ^ - "° \{f + 35^) + {f - 9')vf^ cos ^1 Ef^pf) , ao = ^^^"^''^ (1 + Af_) (9) 



Vdcos^j 
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vr 



where the vector and axial- vector coupling constants are / = 1, 5 = 1.26 and ^fnner — 0-024 represents 
the inner radiative corrections. Integrating over angles one obtains the standard result for the total 
cross section, which can be also related to the neutron lifetime r^, 

= -0 {P + 3,^) Ef)pf^ = 0.0952 (f^V ^ ^ ^o_42 ^ 2^ ^(o)^(o) ^ (^q) 

where /^^^ = 1.7152 is the phase space factor, including the Coulomb, weak magnetism, recoil, and 
outer radiative corrections. For supernova Ve terms of order 1/Mp should be included. The expressions 
for the cross section to that order, including angular distribution, can be found in Ref. [13U . 

For the 'standard' SN with Tp^ = 5 MeV one expects ~8300 e"*" events in Superkamiokande, and 
~360 events in the light water part of SNO. SNO will be able to detect also by the charged 
current reaction on deuterons; one expects about ~80 events of this kind with two neutrons in the 
final state. In KamLAND, which is a scintillation detector, the correlation between the positron and 
the neutron capture 7-rays can be used; one expects ~330 events there. Altogether, it should be 
possible to measure with good accuracy the luminosity and energy distribution of supernova Pe- In 
Superkamiokande the statistics ought to be sufficient to determine also the time dependence of the z^g 
luminosity and temperature. 

It is more difficult to detect z^g since they interact only with neutrons and are expected to have 
lower temperature ( T^^ = 3.5 MeV). Both ^^O (in water Cerenkov detectors) and ^^C (in scintillation 
detectors) have high thresholds for the fg induced charged current reactions, 15.42 MeV and 17.34 
MeV, respectively. Thus, one expects negligible yields for the charged current reactions on these 
targets as long as T^^ is indeed only 3.5 MeV. 

In SNO the 'solar' reaction L'e + d ^ e~ +p + p with mere 1.44 MeV threshold should yield about 
80 events, perhaps sufficient to determine, at least crudely, the temperature T^^ and the corresponding 
luminosity. 

Since cross sections for the charged current Ve interaction with nuclei typically increase quickly with 
Ely, the count rates would increase dramatically if <-> I'x mixing occurs, which is likely to happen. 
Hence observation of the f g signal represents a sensitive test for oscillations. In KamLAND one expects 
only ~2 events for the Ve^'^C — >^^Ngse~ reaction if Tj^^ = 3.5 MeV. That rate increases to ^ 15 for 
vacuum oscillations, and to 27 for the resonant MSW oscillations. (This reaction has an excellent 
signature since one can use the delayed coincidence with the ^^N /3+ decay.) In SuperKamiokande 
the reaction ^^O — >^^F*e~ results in only ~ 20 events for Tj^^ = 3.5 MeV. However, if through 
oscillations the effective T,^^ = 8 MeV, the yield increases dramatically to ~860 events |l^, |l^. The 
electrons from the fg charged current reaction on ^^O can be distinguished, in principle, from the 
positrons from i>e on protons by their angular distribution. 

Lead has been proposed as the target material in OMNIS and LAND supernova neutrino detec- 
tors. The charged current reaction on ^"^Pb induced by has threshold of only 2.9 MeV, but the 
corresponding strength is dominated by the excitation of the giant Gamow- Teller resonance at about 
16 MeV excitation energy. The proposed detectors would register neutrons emitted by the decay of 
the final ^*^^Bi for the charged current reaction or 208p|^* £qj, ^^iq neutral current reaction. Generally, it 
would be difficult to separate the charged and neutral current responses with such scheme. (Although 
with the 'normal' hierarchy, Eqs. (4,5), the neutral current signal would dominate.) However, in 



Ref. [|16[ it was shown that the observation of the double neutrons could serve as a signature of the 
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charged current induced events in the case of oscillations. One drawback is that the corresponding 
cross sections for both reactions are rather uncertain. In fact, the two recent calculations of these 
quantities 17| differ by about a factor of two. Thus, if the lead based supernova detectors are ever 



build, experimental determination of these cross sections will be necessary. 



3 Detecting neutrinos through neutral current scattering 



The supernova I'x, i-e. i^n and Vr with their antiparticles, do not have enough energy to induce 
charged current (CC) interactions. Thus, they can be detected only through their neutral current 
(NC) scattering. In order to detect the NC scattering one has to find, first of all, the appropriate 
signature, i.e. a reaction that can be clearly recognized and separated from the CC channels. Since 
NC scattering is flavor blind, the contribution of the z^e and scattering has to be subtracted in order 
to isolate the effect. This condition more or less eliminates neutrino-electron scattering, where the 
fe and z^e contribution dominates. However, in semileptonic NC scattering the cross section typically 
increases fast with energy, and hence the Vx contribution will dominate the NC yield. (The fact that 
there are four flavors in the I'x flux helps as well.) 

The other difficulty is that in a typical NC reaction there is no spectral information; only the 
number of events per unit time can be measured. Generally, the scattering rate (per s) is: 



dNNc 
dt 

where for SuperKamiokande 



C / dEf{E) 



a{E) 



10 



L{t - At{E)) 



Eb/Q 



C = 9.21 



Eb 



1053 



ergs 



1 MeV 



10 kpc 



-\2 r 



D 



det. mass 



1 kton 



n . 



(11) 



(12) 



T is the spectrum temperature (where we assume {E) = 3.15T, as appropriate for a Fermi-Dirac 
spectrum), f{E), the neutrino energy distribution is in MeV-\ and n is the number of target nuclei 
per water molecule. Also, At{E) is the possible delay caused by the finite neutrino mass. 

Thus, NC scattering rate depends on both the luminosity and temperature, and their effects 
cannot be directly separated. On the other hand, the NC signal is obviously independent of possible 
oscillations between active (as opposed to sterile) neutrinos. 

There are several ways in which the NC signal in existing detectors can be determined. In water 
Cerenkov detectors, one can use the signal proposed in [l^] according to which the Ux neutrinos will 
excite ^^O into the continuum that will deexite dominantly by the emission of either proton or neutron. 
There is a sizable probability (about 30%) that the resulting ^^N or ^^O nucleus will be in a bound 
excited state, as indicated in Fig. |l|. These states, in turn, deexite by 7 emission with characteristic 
energies between 5 and 10 MeV, above the SuperKamiokande threshold, and easily separated from 
the background positrons from i7f,p — > e~^n. In SuperKamiokande one expects about 700 events of this 
kind. 

In SNO the obvious NC signal is the deuteron disintegration, i^x + d n + p + Vx- It can be 
recognized by detecting a single neutron, but no electron (or positron). The rate is obtained by the 
same equations as above, except in Eq. (12) one should replace 9.21 8.28. There will about 400 
NC fx induced events in SNO, with another ~ 85 induced by Ve and Pg- 

In scintillation detectors, such as KamLAND or Borexino, Ux NC scattering with excitation of the 
15.11 MeV T = 1,/'^ = 1+ state in is possible (see Fig. |2|). This process offers a very distinct 
signature and has the further advantage that the corresponding cross section is calculable accurately, 
and has been verified by the KARMEN and LSND experiments. One expects ~ 60 events with 15.11 
MeV 7 in KamLAND. 



5 




15 



N + p 



16 



o 



Figure 1: Schematic illustration of the detection scheme for supernova Vx neutrinos in water Cerenkov 
detectors. 
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Figure 2: Illustration of the excitation of the T = 1, 7=1+ triad in mass A =12 nuclei . 



Finally, as mentioned above, in lead based detectors the single neutron events will be dominated 
by the NC Vx scattering. 

Thus, there will be a rather accurate information on the rate of the NC events. By combining the 
data from different detectors, one can try to determine the Ux luminosity and temperature separately. 

This should be possible, at least crudely, since the mentioned reactions, while all proportional to the 
fx luminosity, will have slightly different dependence on neutrino energy in the various respective cross 
sections. 
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4 Neutrino elastic scattering on protons 



Ideally, one would like to use NC scattering combined with some spectrum information, not just rate as 
in the previous section. As stressed previously also, the seemingly obvious candidate process, neutrino 
- electron scattering, will be dominated by the Vf. and scattering, and thus is not very convenient 
to study the scattering. 

In detectors with low detection threshold, such as the scintillator based KamLAND and Borexino, 
one can, in principle use for this purpose the elastic scattering on protons |^. The corresponding 
differential cross section is 



dcr _ G\Mp 
dTp vr 



(13) 



where cy = 1/2 — 2sin^ 9w = 0.0375, ca = 1.26/2, cm — —fJ-nf^, and ± refers to and u, respectively. 
(We have neglected the possible effect of the strangeness component of the proton). The total cross 
section is proportional to E^, so the signal will be dominated by Ux, particularly above reasonable 
detection thresholds. However, while the recoiling protons scintillate, the scintillation light is quenched, 
compared to electrons or 7. Thus, the relevant observable energies are < 1 MeV, and difficult to detect 
and separate from backgrounds. However, in a sensitive low background detector one might be able 
not only to count the number of events, but actually observe the proton recoil spectrum. The cross 
section, without account of quenching, and for the Fermi-Dirac spectrum of incoming neutrinos, is 
shown in Fig. H. Note the sharp dependence on the neutrino temperatures. 




proton recoil energy T (MeV) 



Figure 3: Cross section of the elastic neutrino scattering on protons for the indicated temperatures of 
the incoming neutrinos. Proton recoil energies without quenching are used. 

^The content of this section is based on the suggestion of John Beacom, for details see [Eil. 
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Let us assume that one will be able to extract from measurement some spectral information on 
the recoiling protons. Would that make it possible to distinguish the cases in which the luminosity 
and temperature conspire in such a way that they lead to the same total number of events, and 
therefore are indistinguishable based only on Eq. (11)? The answer is yes, and how this could 
be accomplished is illustrated in Fig. ^. One can see that the proton recoil spectra sensitively 
depend on the neutrino temperature, with the ratio of the low and high energy yields decreasing 
with the increasing temperature. In a detailed simulation [19j the power of such discrimination was 
demonstrated by taking into account the statistical fluctuation of the expected data. As shown in 
Fig. ^ one expects about 10% resolution on both the temperature and total energy carried by these 
neutrinos. 

It should be stressed once more that the considered NC signal is independent of neutrino oscillations 
into 'active' flavors, i.e. ^ f^,r and obviously ^ v^-- If this signal can be in fact detected, it 
would measure the luminosity and temperature of the hottest component of the supernova neutrino 
emission spectrum. 




0.2 - 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 

Effective energy (MeV) 



Figure 4: Ratio of proton yields, as a function of the effective quenched energy, to the standard case 
of T = 8 Mev, and the total energy emitted in Vx equal to 2x10^^ erg. All considered cases result in 
the same total number of events above the threshold of 200 keV of the effective energy. 
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Figure 5: Monte Carlo simulation of the combined fit to Tj^^ and the total energy carried by such 
neutrinos, E^^. 

5 Conclusions 

In this lecture I have shown how, through the combination of the existing (or soon to be operational) 
detectors, one can determine simultaneously and independently the luminosities and average energies 
(or temperatures) of the three expected components, Ug, and v^, of the next Galactic supernova 
neutrino flux. For a 'standard' supernova near the center of our galaxy, at 10 kpc, I have estimated 
the corresponding count rates, neglecting for a moment the possible effects of neutrino oscillations. 

Having this set of quantities will make it possible to verify, or find deviations, from the basic 
assumptions about the supernova neutrino emission: the equal luminosity in each of the six neutrino 
flavors, and the hierarchy of decoupling temperatures. Also, one should be able to determine the 
total emitted energy, essentially the supernova binding energy, and the total neutrino fluence. Such 
observables will, in turn, severely constrain theoretical models of supernova neutrino emission, and 
allow one to deduce conclusions about the possible role of neutrino oscillations. 

Most of the original results reported here were obtained in a highly pleasurable collaboration with 
John Beacom. The work was supported by the US DOE contract DE-FG03-88ER40397. 



9 



References 



[1] K. Hirata et al, Phys. Rev. Lett. 58 (1987) 1490 

[2] R. Bionta et al., Phys. Rev. Lett. 58 (1987) 1494 

[3] D. N. Schramm and J. W. Truran, Phys. Rep. 189 (1990) 89 

[4] S. van den Bcrgh, Phys. Rep. 204 (1991) 385 ; G. A. Tamman, W. Loffler and A. Schroder, 
Astrophys. J. Suppl. 92 (1994) 487 

[5] J. F. Beacom and P. Vogel, Phys. Rev. D 58 (1998) 053010 

[6] J. F. Beacom and P. Vogel, Phys. Rev. D 58 (1998) 093012 

[7] D. E. Groom et al. (Particle Data Group), Eur. Phys. J. 15 (2000) 1 

[8] J. F. Beacom, R. N. Boyd and T. Mezzacappa, Phys. Rev. Lett. 85 (2000) 3568 

[9] J. F. Beacom and P. Vogel, Phys. Rev. D 60 (1999) 033007 

[10] H. A. Bethe, Rev. Mod. Phys. 62 (1990) 801 

[11] S. E. Woosley et al, Astrophys. J 433 (1994) 229 

[12] H. -T. Janka in Valcano Workshop 1992, F. Giovannelh and G. Mannocchi editors, (Italian 

Physical Soc, Vulcano, 1993) 

[13] P. Vogel and J. F. Beacom, Phys. Rev. D 60 (1999) 053003 

[14] W. C. Haxton, Phys. Rev. D 36 (1987) 2283 

[15] E. Kolbe, K. Langanke, and P. Vogel, to be published 

[16] G. M. Fuller, W. C. Haxton, and G. C. McLaughhn, Phys. Rev. D 59 (1999) 085005 

[17] E. Kolbe and K. Langanke, Phys. Rev. C 63 (2001) 025802 

[18] K. Langanke, P. Vogel, and E. Kolbe, Phys. Rev. Lett. 76 (1996) 2629 

[19] J. F. Beacom, W. M. Farr, and P. Vogel, in preparation. 



10 



